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Recently it has shown that Einstein's field equations can be rewritten into a form of the first 
law of thermodynamics both at event horizon of static spherically symmetric black holes and 
apparent horizon of Friedmann-Robertson- Walker (FRW) universe, which indicates intrinsic 
*Q ! thermodynamic properties of spacetime horizon. In the present paper we deal with the so- 

called f(R) gravity, whose action is a function of the curvature scalar R. In the setup of 
static spherically symmetric black hole spacetime, we find that at the event horizon, the field 
equations of f(R) gravity can be written into a form dE = TdS — PdV + TdS, where T is the 
Hawking temperature and S = Af'(R)/4G is the horizon entropy of the black hole, E is the 
horizon energy of the black hole, P is the radial pressure of matter, V is the volume of black 
hole horizon, and dS can be interpreted as the entropy production term due to nonequilibrium 
thermodynamics of spacetime. In the setup of FRW universe, the field equations can also be cast 
to a similar form, dE = TdS + WdV + TdS, at the apparent horizon, where W = {p — P)/2, E 
is the energy of perfect fluid with energy density p and pressure P inside the apparent horizon. 
Compared to the case of Einstein's general relativity, an additional term dS also appears here. 
The appearance of the additional term is consistent with the argument recently given by Eling 
et al. ( |gr-qc /0602001 ) that the horizon thermodynamics is non-equilibrium one for the f(R) 
gravity. 
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1 Introduction 



The four laws of black hole mechanics within Einstein's theory of general relativity are very 
closely analogous to the four laws of the usual thermodynamics [lj. The quantities that provide 
a basis for the black hole thermodynamics are the Hawking temperature T = k/2it, entropy 
S = A I AG and energy E, where k is the surface gravity, A is the area of the event horizon 
and E = M is the mass of the black hole. The Hawking temperature together with black hole 
entropy is connected through the identity TdS = dE and is usually called the first law of black 
hole thermodynamics [H El [3] . In general, the first law of black hole thermodynamics is related 
with the energy change, dE when a black hole varies from one stationary state to another by 

dE = —dA + work terms, or dE = TdS + work terms. (1) 

The "work terms" are given differently depending upon the type of the black holes. For Kerr- 
Newman black hole family, the first law of black hole thermodynamics states that the differences 
in mass M , the horizon area A, the angular momentum J, and the electric charge Q of two nearby 
black holes are connected by dM = TdS+£ld,J+<&dQ, where f2 = ^j- is the angular velocity and 
$ = Qjq- is the electric potential. The first law of black hole thermodynamics seemingly indicates 
that there should be some relationship between thermodynamics and Einstein field equations 
because black hole solutions are derived from the Einstein equations and the geometric quan- 
tities (horizon area, surface gravity) of the spacetime metric are related with thermodynamical 
quantities (entropy, temperature) of the thermal system. In this regard, Jacobson [3] is the 
first one to explore such an relation. He found that it is indeed possible to derive the Einstein's 
equations from the proportionality of entropy to the horizon area together with the fundamental 
relation 5Q = TdS, assuming the relation holds for all local Rindler causal horizons through 
each spacetime point. Here 5Q and T are the energy flux and Unruh temperature seen by an 
accelerated observer just inside the horizon. For the so-called f(R) gravity, however, Eling, 
Guedens, and Jacobson [5] recently have shown that in order to derive the corresponding field 
equations by using the procedure in [3], a treatment with nonequilibrium thermodynamics of 
spacetime is needed, in which an additional entropy production term dS has to be added to 
the Clausius relation, 5Q/T = dS + dS. More recently, we have shown that a similar entropy 
production term is also required for the scalar-tensor theory, while it is not for the Lovelock 
gravity [6]. 

If one identifies the black hole mass M as an energy E, compared to the standard first law 
of thermodynamics (dE = TdS — PdV), obviously a work term related to the volume change is 
absent in the first law of black hole thermodynamics dE = TdS. More recently, by generalizing 
earlier works |7J, Paranjape, Sarkar and Padmanabhan [8] have considered a special kind of 
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spherically symmetric black hole spacetimes and found that it is indeed possible to interpret 
Einstein's equations as the thermodynamic identity TdS = dE + PdV by considering black hole 
horizon as a boundary of the system, here P is the radial pressure of matter at the horizon. They 
have also shown that the field equations for Gauss-Bonnet gravity and in more general Lancos- 
Lovelock action in spherically symmetric space-time can also be expressed as TdS = dE + PdV . 
These results provide a deeper relation between the thermodynamics of horizon and Einstein's 
field equations of static, spherical symmetric black hole spacetimes. 

The thermodynamical properties of the black holes horizon can be extended to the space- 
time horizons other than black hole horizon. For example, the de Sitter space time with radius 
£, there exists a cosmological event horizon. This horizon, like a black hole horizon, can be 
regarded as a thermodynamical system [9] associated with the Hawking temperature T = l/2ir£ 
and entropy S = Aj AG, where A = Air£ 2 is the cosmological horizon area. For an asymp- 
totic de Sitter space, like Schwarzschild-de Sitter space-time, there still exists the cosmological 
horizon which behaves like a black hole horizon with entropy proportional to the area of the 
cosmological horizon and whose Hawking temperature is given by T = k/2tt, where k is the 
surface gravity of the cosmological horizon. It is easy to verify that the cosmological horizons 
of these space-times satisfy the first law of black hole thermodynamics of the form TdS = —dM 
[lOj . where the minus appears due to the fact that when the black hole mass M increases, the 
cosmological horizon radius decreases. In this frame work, the first law TdS = —dE of cosmo- 
logical horizon thermodynamics does not involve the work term as well. By applying the first 
law of thermodynamics (TdS = —dE) to the apparent horizon of the an (n + l)-dimensional 
Friedman-Robertson- Walker (FRW) universe and by working out the heat flow through the 
apparent horizon, one of the present authors and Kim [TT] are able to derive the Friedmann 
equations of the FRW universe with any spatial curvature, during which one assumes that the 
apparent horizon has temperature and entropy expressed by 

" 1 S-±, (2) 
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respectively, where A is the area of the apparent horizon with radius ta- Also by using the 
entropy expression of a static spherically symmetric black hole in the Gauss-Bonnet gravity and 
in more general Lovelock gravity, they reproduce the corresponding Friedmann equations in each 
gravity. The possible extensions to the scalar-tensor gravity and f(R) gravity theory have been 
studied in reference [12]. In the cosmological setting, related discussions see also, [13 | [T¥ | [15 | [16] . 

On the other hand, we know that the Friedmann equations of FRW universe are the field 
equations with a source of perfect fluid. In this cosmological setup, different from the case 
of black hole spacetimes discussed in [5], there is a well-defined concept of pressure P and 
energy density p. Therefore, it is required to establish a new relationship between the first law 
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of thermodynamics with work term and the Friedmann equations of FRW universe. In this 
respect, the authors of the present work found [17] that by employing Misner-Sharp energy 
relation inside a sphere of radius ta of apparent horizon, it is possible to rewrite the differential 
form of the Friedman equations as a universal form TdS = dE + WdV , where W = (p — P)/2. 
This is nothing but the unified first law of trapping horizon firstly suggested by Hayward [18] . 
We extended this procedure to the Gauss-Bonnet gravity and in more general Lovelock gravity, 
and verified that in both cases, the Friedmann equations near apparent horizon can also be 
interpreted as TdS = dE + WdV, where E is the total matter energy (pV) inside the apparent 
horizon. The thermodynamics of apparent horizon in the brane world scenario also obeys such a 
formula [19]. However, it is interesting to investigate to what extent the horizon thermodynamics 
can be developed in a spacetime against quantum gravitational effects. A feasible way to proceed 
in this direction is to weigh up the effects of higher curvature corrections to Einstein gravity, since 
such corrections naturally occur due to quantum effects [20] . As a special case of considering the 
effect of higher curvature corrections, we will deal with in this paper the so-called f(R) gravity, 
whose action is an arbitrary function of curvature scalar R. When f(R) = R, the Einstein's 
general relativity is recovered. 

This paper is organized as follows. In Sec. 2, we will discuss the field equations of a general 
class of spherical symmetric black hole spacetimes at the horizon and rewrite them into a form of 
the first law of thermodynamics by using Hawking temperature and entropy associated with the 
black hole horizon. In Sec. 3, we will discuss, in the cosmological setup, the corresponding field 
equations for the f(R) gravity. We will show that, at the apparent horizon of FRW universe, 
Friedmann equations can also be written into a form of first law of thermodynamics. But in 
both cases, an additional entropy production term will appear. Finally in Sec. 4 we will discuss 
our results. 

2 Thermodynamic behavior of field equations in the black hole 
spacetime 

We consider a class of static, spherically symmetric spacetimes expressed by the line element 



where U{r) and L(r) are two functions of r. We assume that the function U(r) has a simple 
zero at r = a and U'(a) ^ but has a finite value at r = a. This defines a space-time horizon 
at r = a. Note that L(r) is a finite and continuous function of r, therefore the horizon has a 
finite area A = 4-7rL 2 (a) and hence the associated entropy of the horizon is determined by the 



ds 2 = -U(r)dt 2 + 
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dr 2 + L 2 (r)(d9 2 + sin 2 0# 2 ), 



(3) 
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location of the zero at r = a of the function U(r). The associated Hawking temperature with 
the horizon can be obtained from the periodicity consideration of the Euclidean time. One can 
obtain a Euclidean metric by Wick rotating the time coordinate t —* —it. This metric will be 
regular at the horizon, r = a, if t is taken to be an angular variable with a period (3 = 4ir/U'(a). 
This period is just the inverse Hawking temperature T = 1/(3 of the black hole. So the non-zero 
surface gravity n = U'(a)/2 is related to the Hawking temperature via T = k/2tt. 
The components of the Einstein tensor for the metric ([3]) are given by 

G° = -^(-l + UL /2 + L(U'L' + 2UL")), (4) 
G\ = p(-l + LU'L' + UL' 2 ), (5) 



and 



G\ = G\ = 7 r 7 (2U'L l + LU" + 2UL"), (6) 



where prime stands for the derivative with respect to r. It can be seen readily that the compo- 
nents Gg and G\ of Einstein tensor are not equal but at the spacetime horizon where U(a) = 0, 
they become equal and are given by 

G° \ r =a = G\\ r=a = j^(-l + LU'L'), (7) 

where U' and V are evaluated at r = a. The explicit evaluation of curvature scalar reads 

R = -ttt(-2 + 2UL' 2 + L 2 U" + ALiU'L' + UL")). (8) 

At the horizon, the curvature scalar reduces to 

R\ r =a = --^(-2 + L 2 U" + 4LU'L'). (9) 

Now let us consider the Lagrangian of the f(R) gravity 

L = T^G f{R) + Lm (10) 

where f(R) is a continuous function of curvature scalar R and L m is the Lagrangian density of 
matter fields. Varying the action yields corresponding equations of motion 

f,n(R)R^ - \g^f{R) - VpVufMF) + g^ 2 f, R (R) = ^gt^, (ii) 

where TJJ^ = -^=j^jSt is the energy-momentum tensor for matter fields, V represents the covari- 
ant derivative defined with the Livi-Civita connection of the metric, R^ v is the Ricci tensor and 
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the subscript ", R" denotes the derivative with respect to the curvature scalar R. The above 
equations can be cast to a form 

iV - \R9v = 8nG (j^(R-) T ^ + ^G T ^ r) ' (12) 
where T™ r is a stress-energy tensor for the effective curvature fluid and is given by 

t%T = j^m^iW) - RfMR)) + v^vfAR) - g^ 2 fAR))- (is) 

Setting f, R ( R ) = F ( R ) and usin g the relations V 2 F = ■j=d fJ ,(y/=gg (U 'd 1 ,F) and V^V^F = 
d^d u F — T^ v d\F along with metric (3), one can find the components of T™ r given by 



T° icur) = (\(f(R) - RF(R)) + - j- 2 (2LL>UF> + L W + L 2 UF")j , (14) 

and 

Tl {cur) = -!- (!(/(*) - RF{R)) + (UF" + ^) - -^(2LL'UF' + L 2 C/'F' + L 2 UF")). (15) 

Note that here the prime stands for the derivative with respect to r. From equations (j!4j) and 
(|15jl . it is easy to see that the components j 1 ^ cur ^ anc l j> 1 ( cur ) are no f. e q Ua i m general, but at 
the spacetime horizon, where U(a) vanishes, we have 

T 0(cur) = T l(cur) = -L-^R) - RF(R)) - ^f). (16) 

It is evident from equations (|7|) and (|16|) that the field equations are consistent with f{R) gravity 
at the horizon of the spherically symmetric metric (J3J) provided that the stress-energy tensor of 
matter fields has the form 

2 J(m) =r l(m) j (17) 

at the horizon. 

Note that V 2 F(R) = ^(2LL'UF' + L 2 U'F' + L 2 UF"). At the horizon, namely, U(a) = 0, 
one then has V 2 F(R)\ r=a = U'(a)F' and the — component of equations of motion takes the 
form 

z1± t^ l = 8 - G im T ° m) + ray ^ m - RFm - l u,F ">) • (18) 

which can be rewritten as 

- + ^W- = 4 " (T °° (m) + ^c { \ {m ~ RF{R)) ~ \ u ' f '» l2 - ^ 

Note that here the curvature scalar R is taken its value at the horizon R = R(r)\ r=a . If one has 
a close look at the above equation, it throws light on the notion of horizon entropy, temperature, 
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and energy. So a thermodynamic interpretation of this equation is possible. Multiplying by an 
infinitesimal displacement of horizon da on both sides of this equation, it is trivial to rewrite 
this equation in the form 

_ ^F(R)da + ±U'F(R)LdL = 4^(T ° (m) + j^(f(R) ~ RF(R) - U'F')L 2 da. (20) 

Note that F(R)LdL = ^d(F(R)L 2 ) — ^L 2 F'da. Substituting it into the above equation, we 
finally get 

- ±F{R)da + g d( 4 ^^) ) = TgWfrtfda) + ±(f(R) - RF(R))L 2 da. (21) 

Now we first concentrate on the left hand side of the above equation. We consider the first 
term ^F(R)da = ^f : n(R)da of this equation and convert it to the limit of general relativity 
by setting f(R) = R, we get ^da. This term j^da represents the energy change dE = ^da 
[3 [8] during the infinitesimal displacement da within the Einstein gravity {E = a/2G is just the 
Misner-Sharp energy inside the horizon [T7]; For a Schwarzschild black hole solution, the such 
defined energy E = a/2G gives the Schwarzschild mass). Naturally, we may regard the term 
^jF(R)da = dE as the energy change within the f(R) gravity during a infinitesimal horizon 
displacement da. The second term on the left hand side can be written as ^-d( — x —^^- )- Here 
one can identify ^ = T as the temperature of the black hole and the term ri( 47rL 4( ^^ ) is nothing 
but the entropy change dS of the black hole in the f(R) gravity [2HI22]. So the above equation 
can be rewritten as 

TdS -dE = T° (m) 47rL 2 da + ^(/(-R) " RF)L 2 da (22) 

The first term of the right hand side of above equation is Tq ^ 4irL 2 da, where one can recognize 
47rL 2 as a surface area of horizon with a multiple displacement da of the horizon radius a, which 
corresponds to a change of volume. We write it as ^Kl?da = dV . Since at the horizon one 
has T^ m ' = T^™^ = P, which represents the radial pressure of matter fields at the horizon, 
therefore this term turns out to be the work term PdV against the pressure. Hence the above 
equation can be written as 

TdS -dE = PdV + t Jq( RF )da, (23) 

where A = 47rL 2 (a) is the area of the black hole horizon. It has been shown in references [71 [8] 
that the field equations in Einstein, Gauss-Bonnet and Lovelock gravities can be rewritten as 
the thermodynamical identity TdS = dE + PdV at the black hole horizon of a special class of 
spherical symmetric black hole spacetimes, but here it is evident from the above equation that 
the field equation of f(R) gravity can no longer be cast to the thermodynamical identity TdS = 
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dE + PdV at the black hoke horizon if one defines the horizon energy as E = 1/2G Jq F(R)da. 
It is reminiscent of the argument recently given by Eling et al. . In that reference, they argued 
that in order to derive the filed equations for the f(R) gravity by using relation 5Q = TdS and 
assuming the relation holds for all local Rindler causal horizons through each spacetime point, 
here 5Q and T are the energy flux and Unruh temperature seen by an accelerated observer 
just inside the horizon, and S = f rA/AG is the horizon entropy, one has to add an additional 
entropy term to the Clausius relation so that 5Q/T = dS + dS, here the expression of dS is given 
in [5]. The appearance of the additional entropy term is interpreted as entropy production term 
in the nonequilibrium thermodynamics. That is, in their setting, the horizon thermodynamics 
for the f(R) gravity is a nonequilibrium one. In the spirit of [5], it is then natural to regard the 
last term in (I23p as the entropy production term in our setup. In other words, if one includes 
the effects of higher curvature correction to the Einstein gravity and the curvature scalar R in 
general relativity is replaced by a generic function f(R) then the usual thermal behavior of the 
field equations at black hole horizon [8] will not be sustained by the system but will shift from 
an equilibrium setup to a non-equilibrium one [3]. Instead of satisfying the thermodynamical 
identity TdS = dE + PdV at the black hole horizon, we have to put in the effects of higher 
curvature correction to the horizon entropy grown up internally because of the system being 
out of equilibrium [5] . So one may have a balance thermal identity for higher curvature gravity 
theory 

dE = TdS - PdV + TdS (24) 

where TdS is just the last term of (|23p . which is the entropy production term added to balance 
the inequality TdS > dE + PdV . We note that the entropy production term for the f(R) gravity 
has the form 

as^i^m^ (25) 

in our setup of spherically symmetric black hole horizon. It is evident from this equation that 
the entropy production term vanishes at the limit f{R) = R of the Einstein's gravity, and the 
usual thermal behavior TdS = dE + PdV of the field equation at the horizon is recovered. 

Clearly the above interpretation of nonequilibrium thermodynamics of black hole horizon is 
not unique. One may combine the last term in (I2ip to the first term so that one can define a 
new energy associated with the black hole horizon as 

dE = h Fda + h {f{R) ~ RF ^ L2(la - ( 26 ) 

In this way, we can rewrite (|21|) to the standard first law of thermodynamics 

dE = TdS - PdV (27) 
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Take an example, let us consider a constant curvature black hole solution of the f(R) gravity 
without matter fields. In that case, we have from that f^Ro = 2f(Ro), where the i?o is 
the curvature scalar of the black hole solution. The field equations (llip then can be cast to the 
standard Einstein's equations with an effective cosmological constant A = 3/£ 2 = -Ro/4- Thus 
we have a simple Schwarzschild-de Sitter (anti-de Sitter) black hole solution depending on the 
sign of the effective cosmological constant, which has the metric functions, L(r) = r and 

2GM r 2 , , 

U(r) = l-—--^, (28) 

where M is the mass of the black hole in the Einstein gravity. In terms of the black hole horizon 
radius U(r)\ r=a = 0, the mass can be expressed by 

„2\ 




On the other hand, it is interesting to note that integrating (|26p yields 

E = F(R )M. (30) 



Clearly, the black hole mass M does not satisfy the equation (I27p with P = 0, but E does. 
Therefore defining E as the horizon energy of the black hole is reasonable. 



3 Thermodynamics of apparent horizon of FRW universe 

In the previous section we have discussed the behavior of the field equations for a static, spher- 
ically symmetric black hole spacetimes at the black hole horizon. In this section we generalize 
the thermodynamics of black hole horizon to a dynamical apparent horizon in a FRW universe 
within the f(R) gravity. We consider an (n+l)-dimensional, spatially homogenous and isotropic 
FRW universe described by the metric 

dv 2 

ds 2 = -dt 2 + a{t) 2 {- — 1 +r 2 dti 2 n „ 1 ), (31) 

where a(t) is the scale factor of the universe with t being cosmic time and dfl 2 l _ 1 is the metric 
of (n — l)-dimensional sphere with unit radius and the spatial curvature constant k = 1, and 
— 1 correspond to a closed, flat and open universe, respectively. Using the spherical symmetry, 
the metric (|3T1) can be rewritten as 

ds 2 = h ab dx a dx b + f^O^!, (32) 

where f = a(t)r and x° = t, x 1 = r and the two dimensional metric h a i, = diag(— 1, a 2 /I — kr 2 ). 
The dynamical apparent horizon is determined by the relation h ab d a fd},f = 0, which implies that 
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the vector Vf is null on the apparent horizon surface. The explicit evaluation of the apparent 
horizon gives the apparent horizon radius 

i = # 2 + 4>> (33) 
r A a z 

where H denotes for the Hubble parameter. The associated temperature T = k/2tt with the 
apparent horizon is defined through the surface gravity 

1 



d a (V-hh ab d b r), (34) 



2y/-h 

which gives 



"--^S?- (35> 

From this equation one can see that when < 2HrA, the surface gravity k is negative. If 
one still uses T = k/2tt to define a temperature of the apparent horizon, we obtain a negative 
temperature! This case is quite similar to the case of cosmological horizon for the Schwarzschild- 
de Sitter spacetime. In this case, one should use T = \n\/2n to define the temperature of the 
apparent horizon. 

The total matter energy inside a sphere of radius ta of the apparent horizon is given by 

E = Q n f>, (36) 

where V = 0, n f A is the volume of the sphere within the apparent radius, and p is the energy 
density of the perfect fluid in the universe. We consider FRW universe as a thermodynamical 
system with apparent horizon surface as a boundary of the system. In general the radius of the 
apparent horizon fj\ is not constant but changes with time. Let dr^ be an infinitesimal change 
in radius of the apparent horizon during a time of interval dt. This small displacement df^ in 
the radius of apparent horizon will cause a small change dV in the volume V of the apparent 
horizon. Since the energy E inside the apparent horizon of FRW universe is directly related 
with the apparent radius ?a therefore an infinitesimal change dfj^ will result in a change dE in 
energy E and is given by 

dE = nQ, n r%pdr A - n^l n f n A {p + P)Hdt. (37) 



Here, we have utilized the continuity equation (j39jh 

Now we consider the field equations (|12|) . the stress-energy tensor T™ is taken to be the one 
of the perfect fluid with time dependent energy density pit) and pressure P(t) 

T™ = {p + P)U^U v + P 9flv , (38) 
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where Uu is the four velocity of the fluid. The conservation of the stress-energy tensor, j:^" 1 ) = 
0, leads to the continuity equation 

p + nH{p + P) = 0, (39) 

where the dot represents derivative with respect to cosmic time t. 

Using the relations V 2 F = -^d^^gg^ d v F) and V^V U F = d^d u F-T^ u d x F along with 
the (n + l)-dimensional metric (|31|) of the FRW universe, one can find the components of T™ r 

T ™ = -p^(l(f(R)- RF (R)) + nHF R (R)R), (40) 
T^ r = — !— Q(/(#) ~ RF(R)) + F R (R)R + F AR (R)R 2 + (n - l)HF R (R)R^j g ti (41) 
where F(R) = df /dR, F iR = dF/dR and F iRR = d 2 F/dR 2 . The field equations become 

k " ? = -<^W(R) ( ( " + p) + ^a ( d{F * {R)R) - H} "^ R )) • (43> 

where d(F tR (R)) = d{F )R {R)) / dt. These are two Friedmann equations for the f{R) gravity |23j . 
Note that, as mentioned in the previous section, in this gravity theory the black hole entropy 
has a relation to its horizon area S = -^F(R) (Horizon- We assume it also holds for the apparent 
horizon. Using relation H 2 + \ = Jr and taking differential of it, one gets 

a r A 

- \d? A = H(H - ^)dt. (44) 



r A a 



Substituting it into (|43l) and then multiplying in both sides with a factor nQ n r A , we reach 

-V(n(n - l)n n rl' 2 F(R)df A )/AG = nSl n r n A H ( (p + P) + J-(d(F R (R)R) - HF R R)) dt. 

(45) 

Now consider the differential 

dinting 1 F(R)/ AG) = n(n - l)Q n r A - 2 F(R)dr A /AG + n^l n f n ^ 1 F fi {R)Rdt/AG, (46) 
and substitute it into (|45|) . one gets 

1 _ d{ na n r n A 1 F{R) ) = A + p) + 1 (d(F(J2)A) _ HF iR R)) dt 



2vrf A 4G " 71 V 8vrG 



nQ n r n A - 2 F R Rdt 

+ — t^ 5 — (47 > 
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Multiplying both sides of this equation by a factor —(1 — f A /2HrA) and identifying that the 
surface gravity n = (— l/f^)(l — r A /2Hf : A ) at the apparent horizon, one can obtain 



^ d ^r) = ^-^)n^nflH^p + P) + ^(d(F R (R)R)-HF R Ryt 

-dt, (48) 



r _ r_A ntt n ^- 2 F R R 
2Hf A > 8vrG 



where A = nQnf 1 ^ 1 is the area of the apparent horizon and VL n = ir n / 2 /T(n/2 + 1) being the 
volume of an n-dimensional unit ball. In addition, let us note that all quantities in the above 
equation are evaluated at the apparent horizon. Recognizing k/2it as the temperature T of the 
apparent horizon and the quantity inside the parentheses on the left hand side of this equation, 
is just the entropy S = ^F(R)\r A °f the apparent horizon in the f(R) gravity. Hence, having 
identified temperature T and entropy S , the above equation can be written as 

TdS = -nn n f 7 X(p + P)Hdt + ^n n f 1 \- 1 (p + P)dfA 

+T A (Hf\{d{F tR (R)R) - HF R (R)R) + F iR (R)R) dt. (49) 

Now we turn to the total matter energy E = Q, n f\p inside the apparent horizon. Using (137)) . 
the above equation becomes 

TdS -dE = - X -{p- P)d^l n f n A )+T^ (Hr 2 A (d(F^ R {R)R) - HF R (R)R) + F r {R)R) dt. (50) 

Note that {p — P)/2 = W is the work density |18j . Hence the equation (|50p can be further 
rewritten as 

dE = TdS + WdV - T^- (Hr 2 A (d(F tR (R)R) - HF, R (R)R) + F R (R)R) dt. (51) 

In the cosmological setting, we have already shown [17] that the field equations obey the universal 
form dE = TdS + WdV at the apparent horizon of FRW universe in Einstein, Gauss-Bonnet and 
Lovelock gravities. But it no longer holds for the scalar-tensor theory [6]. Here we see from the 
equation (|51|) that the universal form dE = TdS + WdV also does not hold for the f{R) gravity 
at the apparent horizon of FRW universe because of the appearance of the additional term 
T^ (Hf\(d(F R {R)R) - HF R (R)R) + F R (R)R) dt. Of course when f(R) = R, the additional 
term is absent, the usual result for the Einstein gravity is recovered [17 . In the spirit of the 
argument of [5], once again, the additional term can be interpreted as follows. This additional 
term is developed internally in the system as a result of being out of equilibrium. In order to 
get a balance relation for a universal form, one may take 

dE = TdS + WdV + TdS (52) 
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The term dS is the entropy production term grown up internally due to the non-equilibrium 
setup within the f(R) gravity at the apparent horizon of FRW universe. Comparing equation 
dSH) and ([52]), we get 

dS = ~ (Hf\(d(F R (R)R) - HF R (R)R) + F R (R)R) dt. (53) 

If one further defines dS = d{S + S) as an effective entropy change during the infinitesimal 
displacement dr^ of the apparent horizon of FRW universe in interval dt, the above equation 
(|52|) can be rewritten as 

dE = TdS + WdV, (54) 
where S = S + S is the effective entropy associated with the apparent horizon of FRW universe. 

4 Conclusion and Discussion 

Black hole thermodynamics implies that there must exist some relation between gravity theory 
and the laws of thermodynamics. Indeed, it has been shown that at the horizon of spherical 
symmetric black hole spacetime, field equations can be written into the form of the first law 
of thermodynamics, dE = TdS — PdV, not only for the Einstein gravity, but also Gauss- 
Bonnet gravity and more general Lovelock gravity 0(8]. In the cosmological setting, applying 
the Clausius relation 5Q = TdS to the apparent horizon of a FRW universe, one can derive 
corresponding Friedmann equations of the universe with any spatial curvature in the Einstein, 
Gauss-Bonnet and Lovelock gravity theories Further we have shown that at the apparent 
horizon, the field equations of gravity can also be cast to a similar form of the first law, dE = 
TdS + WdV , in the Einstein, Gauss-Bonnet, and Lovelock gravity theories [17]. To what extent 
does such a formulism hold? In the present paper we have discussed the so-called f(R) gravity. 

In the setup of static, spherically symmetric black hole spacetime, we have shown that the 
field equations of the f(R) gravity can be rewritten as dE = TdS — PdV + TdS, at the black 
hole horizon. Here E = 1/2G Jq F(R)da is regarded as the horizon energy; when f(R) = R, the 
horizon energy reduces to the one given in [7J for Einstein gravity. T is the Hawking temperature 
of the black hole and S = F{R)A/AG is the horizon entropy of the black hole, and P is the 
radial pressure of matter fields at the horizon. The additional term dS given in (j25[) can be 
regarded as an entropy production term due to a nonequilibrium thermodynamics setup [5], 
where it is argued that in the setup of a class of Rindler causal horizons, in order to derive 
the field equations for the f(R) gravity, an entropy production term has to be added to the 
Clausius relation 5Q = TdS + TdiS. However, in our case, we have argued that the additional 
term can be combined to the horizon energy term by defining a new horizon energy (I26p . In 
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this way, the field equations of f(R) gravity can be cast into the standard form of the first 
law of thermodynamics (|27|) : dE = TdS — PdV, which reduces to the first law of black hole 
thermodynamics when P = 0. 

In the setup of cosmology of FRW universe, compared to the cases of Einstein gravity, 
Gauss-Bonnet gravity and Lovelock gravity, we have found that an additional term appears 
when the Friedmann equations of the f(R) gravity are rewritten into a form of the first law of 
thermodynamics at the apparent horizon of the universe, dE = TdS + WdV + TdS. As usual, 
here E is the total energy of matter inside the apparent horizon, T and S are the associated 
temperature and entropy with the horizon, and W = {p — P)/2 is the work density. The 
additional term in this case is given by (|53j) . Once again, the additional term can be interpreted 
as the entropy production term associated with the apparent horizon in the nonequilibrium 
thermodynamics within the f(R) gravity in the spirit of [5]. 

The authors of [5] used the modified Clausius relation 



to a Rindler horizon of spacetime with the assumption that the horizon has a temperature 
T = 1/2-7T and entropy S = AF{R)/AG, where A is the horizon area. They defined the heat as 
the mean flux of the boost energy current of matter across the horizon 



In this frame work, in order to get the correct field equations of the f(R) gravity, they worked 
out that the entropy production term diS is 



where 9 = d(lnd 2 A)/d\ is the expansion of the congruence of null geodesies generating the 
horizon and A is an afflne parameter. 

We note from (f25j) . (f53l) and (f57|) that in the different settings, the additional entropy term has 
a different from for the f{R) gravity. The same happens in the scalar-tensor gravity [BJ. It would 
be of great interest to further understand this issue. In addition, why does such an additional 
term appears for the f{R) gravity and scalar-tensor gravity, while it does not in the Einstein's 
gravity, Gauss-Bonnet gravity and more general Lovelock gravity. This is a quite interesting 
question, which deserves further investigation. This might be related to the observation that 
the field equations for the Einstein's gravity, Gauss-Bonnet gravity and Lovelock gravity can 
also be derived from a holographic surface term [24J. 



5Q = TdS + TdiS, 



(55) 




(56) 




(57) 
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